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Abstract 



We review the currrent cosmic parameter determinations of relevance to inflation using the 
WMAP-lyear, Boomerang, CBI, Acbar and other CMB data. The basic steps in the pipelines 
which determine the bandpowers from the raw data from which these estimations are made are 
^ ' summarized. We forecast how the precision is likely to improve with more years of WMAP 

Jf[j . in combination with future ground-based experiments and with Planck. We address whether 

the current data indicates strong breaking from uniform acceleration through the relatively 
V^ ' small region of the inflaton potential that the CMB probes, manifest in the much-discussed 

^P , running spectral index or in even more radical braking/breaking scenarios. Although some 

weak "anomalies" appear in the current data, the statistical case is not there. However increased 
■^ ' precision, at the high niultipole end and with polarization measurements, will significantly curtail 

f^ , current freedom. 
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to appear in: Int. J. Theor. Phys. 2004, ed. E. Verdaguer, "Peyresq Physics 8", "The Early 
Universe: Confronting theory with observations" (June 21-27, 2003) 

I. Introduction to CMB Power Spectra 

A Overview 

The three Peyresq lectures covered CMB theory and analysis. The main content and most of the 
relevant references are given in Bond, Contaldi and Pogosyan (2003, hereafter BCP). Although a 
summary of that material will be given here, in this paper we will emphasize the impact on inflation 
phenomenology of the experiments WMAP, Boomerang, CBI and Acbar, in conjunction with DASI, 
VSA, Archeops, Maxima, TOCO, earlier CMB experiments and with large scale structure (LSS) 
observations. We will also remind the reader of the great success now of CMB determinations of 
the material content of the universe, including dark energy and dark matter. A detailed treatment 
of LSS and the relation to the Sunyaev-Zeldovich (SZ) effect, the cluster-dominated upscattering of 
CMB photons from hot electrons in the cosmic web, was covered in the third lecture, and in Bond 
et al. (2004), Readhead et al. (2004); see also Bond and Crittenden (1999). At Peyresq, Page 
looked to the future with WMAP and the Atacama Cosmology Telescope (ACT) and de Bernardis 
covered Boomerang and beyond, concentrating on polarization. We will also explore here how such 
future CMB observations may help to further discriminate among inflation models. 

The theory of Gaussian primary anisotropics, those arising from linear physics operating in the 
early Universe, is in good shape. For the current data, speedy codes efficiently using past-history 
integrations such as CMBfast and CAMB are quite adequate, and have been "validated" with 
codes solving hierarchies of multipole equations. As precision in the data improves, fresh looks at 
all aspects of the accuracy are worthwhile, and are being done. 

For secondary anisotropics, those arising once nonlinearity develops, the computational state of 
the art continues to need much further effort. This includes the important component which rears it 
head at small angular scales associated with the SZ effect. The statistically inhomogeneous Galactic 
foregrounds offer even more of a theoretical challenge, and this has only partly been addressed. The 
direct interface with observations for these many non-Gaussian signals is much more complex than 
for Gaussian primary anisotropics. Because all the signals are superimposed, the separation of the 
primary, secondary, foreground and extragalactic components inevitably complicates the move from 
multifrequency CMB data to determination of cosmic parameters from "cleaned" primary CMB 
power spectra. 

Major efforts by many groups around the world have been put into developing the statistical 
pipelines which process and clean the raw CMB data to allow efficient and accurate confrontation 
with the theory. Even if primary power spectra are the primary target to be extracted from the 
data, the case so far, the sophistication level required is high: processing timestreams or interfer- 
ometer visibilities, making maps in position or "momentum" space, filtering, cleaning, separating 
component signals, compressing, always with new tools to explore systematic effects and anomalies 
that inevitably appear. The step from raw data to primary power spectrum is enormous, from 
power spectrum to parameters small. Most of the developments were driven by the compelling 
necessity of the CMB teams, consisting of theorists/analysts as well as experimentalists, to extract 
accurate science from beautiful large datasets such as Boomerang and CBI, yet remain within com- 



putational feasibility. Polarization analysis is receiving much further attention now, and in spite of 
the many algorithmic advances in the analysis pipeline in recent years, many more are needed for 
the rosy forecasts of high precision described here to be realized. 

B Grand Unified Spectra of CMB data 

Fig. n shows a "Grand Unified Spectrum" Ci as a function of multipole i compressed on to a small 
number of bandpowers, derived from a June 2003 compilation of CMB data. A best-fit inflation- 
based uniform-acceleration model to that data is also shown. In place of determining "cosmological 
parameters" from the observed bandpowers for various experiments, the parameters here are in 
fact power amplitudes in multipole bands chosen by the analyst. Other parameters characterizing 
experimental uncertainties in calibration and beams are also determined simultaneously. In deter- 
mining errors on the cosmological parameters and/or the GUS bandpowers, the other variables are 
marginalized (probability distributions are integrated). 

To see the remarkable evolution that has occurred in Ci in just a few years we refer the reader to 
Fig. 1 of BCP. This shows a sequence of 4 GUS derived for the data extant in JanOO, Jan02, Jun02 
and Jan03. (For experimental acronyms and much more detailed discussion of how the results from 
the different experiments were used see BCP.) The GUS are in excellent agreement with each other 
and with the first-year WMAP results, hereafter WMAPl. All pre-WMAPl CMB results rehed 
heavily on COBE's DMR results, which anchored the low i < 30 multipoles. Fortunately WMAPl 
spectacularly confirmed DMR. JanOO included TOCO and the Boomerang North America test 
flight, as well as 19 previous experiments, including upper limits. Jan02 included the Boomerang 
long duration balloon flight, DASI and Maxima as well. Jun02 had CBI and VSA added as well. 
By JanOS, improved Boomerang results were added to preliminary 2-year CBI results, extended 
VSA results, Acbar and Archeops. 

By Mar03, the WMAPl data were incorporated, and recalibrated CBI 2 year and VSA results 
were included. The recalibration was tiny but the errors on the calibration decreased by a factor 
of 2. The latter was done off WMAPl, using observations of Jupiter. In GUS methods the 
relative calibrations and their errors come out as a byproduct. BCP showed these were in excellent 
agreement with cross-comparisons made by Eric Hivon between WMAPl and Boomerang maps, and 
with the CBI calibration using WMAPl. A JunOS compilation utilized GUS-based recalibrations 
of Boomerang and of ACBAR. The JunOS GUS shown in Fig.^differs only slightly from the MarOS 
GUS used in BCP. This is basically because experimental calibrations and beam sizes are internally 
determined and marginalized over in making the GUS: if the method is correct it will give robust 
results - and it does. The CBI 2 year data has now been released (Readhead et al. 2004). The 
JunOS GUS in Fig. ^ does not include the new VSA data (Dickinson et al. 2004), but with all the 
other data it does not change the spectra or cosmological parameter determinations much. 

WMAPl dominates the i < 600 bands. Unless explicit joint analyses are done of experiments 
with significant overlap with WMAP's all-sky coverage, the bandpowers for such experiments should 
be dropped from the GUS and parameter determinations. Thus COBE and Archeops were not 
included in the MarOS and JunOS compilations. 
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Figure 1: An optimal Grand Unified Spectrum for tlie post-WMAP Jun03 data is shown 
as crosses. Tliis GUS is a maximum-likeliliood determination of the power in 26 (top-hat) 
bands, with calibration and beam uncertainties of the various experiments fully taken into ac- 
count. A best-fit inflation-based uniform-acceleration ACDM C^ spectrum to this Jun03 data 
is shown as solid green. The parameters are {Otot, ^A, f^fo/i^, ficdm^^, "-s, Tct^^h^a^} = 
{1.0,0.712,0.0224,0.118,0.957,0.108,13.7,0.698,0.84}. The WMAPl data optimally compressed 
on to 49 bandpowers is also shown for low i as triangles in the lower I part of the figure, to high- 
light the two low-£ "anomalies" , at ^ of a few and at £ ~ 20. The best-fit C^ would fit better with 
a slight downward tilt below ^ ~ 30 and beyond £ ~ 500, which a scale-dependent ns{k) could do 
(§nj). Peak and dip locations derived from this optimal spectrum are shown in Fig. [3 



C Peaks, Dips and Damping 

Fig. 121 illustrates some of the "pillars" that we were looking for in the TT data to verify the 
paradigm, exemplified by the best-fit power-law model. (1) The effects of a large scale gravitational 
potential at low multipoles, manifested by a "Sachs- Wolfe" plateau; note the integrated SW upturn 
associated with JIa driving a time-varying gravitational potential at low i and the upturn at higher 
i associated with photon compressions and rarefactions. 

At higher i the next pillars are (2) the pattern of acoustic peaks and dips and (3) the damping 
tail. These are governed by the comoving sound speed r^ = 146 it 3 Mpc and damping length 
Rd = 10.3 lb 0.21 Mpc at photon decoupling, and are scaled with the angular-diameter distance 
relation TZdec to decoupling to get the associated ig and £/) values which define the peak/dip and 
damping structure. 

Fig.EJshows the peak/dip locations and amplitudes and their 1-sigma error bars, as determined 
in BCP, contrasting Jan03 with Jun03. The exercise in BCP was to do this directly from the 
relatively broad-band GUS, using a model that slides a sequence of bands across the data sets. 
In spite of the coarseness of the bands, the peak/dip results are highly accurate as long as all 
band-to-band correlations are taken into account. These peak and dip parameters have also been 
determined for individual experiments such as TOCO, Boomerang, CBI, Archeops, and of course 
for WMAPl, the latter described in Page et al. (2003). Values are given in BCP, Table 2. 

There is good evidence from WMAPl, Boomerang and other CMB analyses that the statistics 
of the primary anisotropics are predominantly (4) Gaussian, i.e., have maximal randomness for a 
given power spectrum. Finding (5) secondary anisotropics associated with nonlinear phenomena, 
due to the SZ thermal and kinetic effects, inhomogeneous reionization, weak lensing, etc. is expected. 
The anomalous extra power at high I over the best-fit primary model evident in Fig. ^ arises from 
combining CBI and Acbar data. Assuming an SZ spectrum makes the case somewhat stronger, 
and suggests this pillar may have been seen (see Fig. Eland Bond et al. 2004, Goldstein et al. 2003. 
Readhead et al. 2004). 

(6) Polarization must be there, with forecasted C\ -patterns of peaks and dips intimately 
related to, though with different phases than, those for TT. As well there is a specific peak/dip 
pattern in the TE cross-correlation of the E-mode with the total intensity predicted. The current 
status of polarization is a broad-band EE detection consistent with inflation by DASI and of course 
the remarkable TE cross correlation of WMAPl. More EE detections are expected soon from 
WMAP2, CBI and Boomerang, and there are many planned experiments that should exquisitely 
determine the EE and TE spectra. Figs. |31 lU and El show the TE and EE spectra for the best-fit 
model and forecasts of polarization detections by WMAP4, by Planckl, and by a fiducial ground- 
based large-array polarization-sensitive telescope. 

A future goal for CMB researchers is to find (7) the anisotropics induced by gravity wave 
quantum noise. Not all inflation models predict this, but the well known simple relation between 
tensor tilt nj, the deceleration parameter q and the tensor-to-scalar power ratio shown in Fig. [21 
suggests there may be a strong enough signal to detect. For the Jun03 best-fit model shown with 
scalar index Ug = 0.957, a uniform acceleration model yields a tensor-to-scalar contribution of 0.17 
c/. the Spergel et al. (2003) upper limit of 0.36 at the 95% CL. A holy grail for the subject is to 
detect the B-mode of polarization at low L For these best-fit parameters, it may be do-able with 
the Planck satellite as the lower panel of Fig. JH illustrates. A nice figure summarizing EE and BB 
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Figure 2: Various pillars which determine the inflation-based uniform-acceleration Ci spectrum 
shown that best-fits the Jun03 data are highlighted. The low i part of the figure repeats as crosses 
the two Jun03 bands of Fig. ^ and contrasts these with the equivalent two bands for the pre- 
WMAPl Jan03 data, which includes COBE and Archeops. The low i shape arises both from the 
climb of the photons through the potential well at decoupling and the propagation of the photons 
through a time varying gravitational potential along the line of sight. Also shown at the bottom 
is the expected tensor component for uniform acceleration with tensor tilt nt = Ug — I = —0.043, 
and above it the Spergel et al. (2003) 95% CL limit, corresponding to Vqw IT^C, < 0.36. The higher 
£ part of the figure shows the peak/dip locations (pk/dipj ^^^ heights Cp^./^ipj, as determined 
from the BCP maximum likelihood sliding-band approach. The points with slightly larger errors 
are for the pre-WMAPl Jan03 data and the heavier are for the post-WMAPl Jun03 data. The 
triangles show the values obtained by ensemble-averaging peak/dip locations and heights over the 
large C^-database used in BCP. (Only a weak prior was applied, which allows large movement of 
peak locations associated with the geometry, hence is preferable to more restrictive priors for this 
application. Note how well these statistically-averaged peaks and dips match those of the best-fit 
model.) 
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Figure 3: Idealized forecasts for detections of (top to bottom) TT, TE, EE and BB bandpowers 
by WMAP with four years of 94 GHz channel data. This is a conservative estimate given the 
other WMAP frequencies, but idealized in the sense that foregrounds and other experimental 
complications are ignored. Simulations use the best-fit Jun03 model shown. The JunOS GUS are 
the crosses in the TT panel. Note the very different scales on the low and high d. panels of the figure. 
As well as the target scalar C^, the tensor C\ contributing at the level predicted if n^ = n^ — 1 is 
also shown. The tensor shape is repeated with the amplitude corresponding to the current WMAPl 
upper limit (dashed) on the tensor to scalar ratio. Primary scalar perturbations have no B-mode 
in linear perturbation theory, but they are induced by lensing (Hu 2000). These and the unknown 
foregrounds present severe challenges for confirmation of the gravity-wave induced B-mode. The 
B-mode bandpower spacings differ from the spacings used for the upper panels. 
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Figure 4: Forecast for how well Planck can do with just its 150 GHz channels for one year of data, 
quite a conservative estimate. Even so, note the anticipated detection level of EE and TE at low 
I, sharpening the tc determination, and the possibility of a statistically significant direct tensor- 
induced B-mode detection. Although ground-based experiments at high resolution should have a 
huge pre-Planck impact (e.g.. Fig. |SJ, the all-sky nature of Planck, its large set of polarization- 
sensitive frequencies and likely longer observing time than that assumed here will make it extremely 
powerful to sort out the many signals that complicate the "primary" quest. 
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Figure 5: Forecast for how well a high-resolution ground-based polarization-sensitive large-array 
experiment like ACT or SPT can do, in conjunction with WMAP4. The specific assumptions used 
are given in the caption of Fig. ^J and are conservative over what might be achievable from the 
ground using bolometer arrays. There is a trade off of sky coverage and noise, {fsky = 0.024 
here, with larger and shallower often improving cosmic parameter determinations, and smaller and 
deeper making the lens-induced B-modes potentially detectable.) Apart from the current DASI, 
WMAP, Boomerang and CBI data already in to flesh out the EE spectrum, a number of other 
experiments are planned. These include some in the very near future, e.g., QUaD and BICEP. 
Forecasts for the proposed QUIET with HEMT arrays look as promising as those for polarized 
ACT/SPT. Given the expected signal levels, all ground and space CMB information available will 
be needed and used to get clean primary polarization results. 



bandpower forecasts for various experiments such as QUaD and BICEP in comparison with Planck 
is given in Hivon &: Kamionkowski (2002). Our forecasts for these are in accord. 

The very tiny B-mode signal predicted and the unknown nature of the polarized foregrounds 
would make this a very hard task indeed, one happily defining a long term future for CMB research 
as the community plans a future NASA CMBPol satellite as the next step in space after Planck. 

II. Cosmic Parameters &; Inflation Issues 

A The Conventional Minimal Parameters 

To the well known minimal inflation-motivated set, {i}f)]i'^,i^cdm^'^j^ki^A,ns,Tc, Ag}, defining 
the allowed Ci shapes, other parameters are sequentially added to probe more complex models 
of inflation or of matter content. {0,j is the density parameter pj relative to the critical density 
Per = 10.5 h^ kev cm~'^, where h is the Hubble parameter in units of 100 km s^^ Mpc~ , so ilfeh^ 
and r^cdmh^ are proportional to the physically more relevant baryon and cold dark matter densities.) 
The vacuum energy is Q\ and the mean curvature energy for our Hubble patch is il^, in terms 
of which the total energy content is i^tot = 1 — f^fc- Sample content extras include a subdominant 
light neutrino component, enhanced relativistic particles ilgrh^, e.g., as products of decay of massive 
particles, adding dynamics to f^A, e.g., through a parameterized time-dependent pressure-to-density 
ratio, WAit). 

Astrophysical complications associated with the "late time" impact of reionization of the Uni- 
verse, presumably by the injection of energy from the first stars, are encoded in a single parameter 
Tc, the depth to Compton scattering when the Universe reionized. tc could have a complex tem- 
poral and spatial structure, and, although the CMB primary anisotropics are not that sensitive to 
the details, finding such signatures in the data is a target of planned high resolution experiments. 

B The Inflation Parameters 

Inflation fluctuations are assumed to have a Gaussian statistical distribution, fully specified by a 
power spectrum of curvature perturbations "Pi^^om (^)- The minimal set has only two parameters, the 
overall initial power spectrum amplitude Ag = V^^^^ikn), evaluated at a normalization wavenumber 
kn, and a single spectral scalar index ns{k) = ns{kn). 

There is of course a vast literature on perturbation theory as applied to inflation (e.g.. Bond 
1994, 1996 for the approach described here). Basic variables are the inflaton field Scpinf', other scalar 
field degrees of freedom 5(j)is which can induce scalar isocurvature perturbations; two gravitational 
wave polarization modes h+,hx- One can encode scalar metric perturbations and their variations 
through the inhomogeneous scale factor a(x, t), Hubble parameter H{x, i) and deceleration param- 
eter q{x,t) = —dlnHa/dlna. For example, the scalar 3-curvature is — 4V^lna. Infiation ends 
when q passes from negative to positive. 

Certain time hypersurfaces upon which to measure the perturbations simplify things quite a 
bit. Sample choices are 4>inf, Ina, InH, In(Ha) and conformal time r = J d\na/{Ha). The 
power spectrum Vinaln (^) ^^ scale factor fluctuations evaluated on time hypersurfaces on which 
the Hubble parameter H is constant becomes time-invariant for wavenumbers k/Ha <C 1 (outside 
of the "horizon"), if there is just one dynamically important scalar held, and remains so until 
fluctuations regain causal contact. If the universe has no net mean curvature, (fcom = (^Inal//, , 



measuring the curvature on comoving hypersurfaces. Another variable used extensively is (, which 
reduces to ipcom if k/Ha <^ 1, hence we sometimes refer to V( in the figures. 

The gravity wave power Vcwik) = {k^ /2tt'^) < h\ + h\ > = {k^ /2tt'^) < hijh^^ > /2 used 
here and in Bond (1994, 1996) is defined to be consistent with the conventions of GW detection 
research, with GW mode functions being the usual hx = hi2 = /121, h+ = {hn — /i22)/2. Most 
people in inflation define a Vh{k) = 2VGw{k), which is what At = Vhikn) is defined in terms of. 

Quantum fluctuations in gravity waves must occur during inflation. The only question is how 
T^Gwik) compares with Vip^^^{k). Thus the minimal Ag, Us should at the very least be augmented 
by an At and rit- As well subdominant isocurvature components can be added, Ais and nj^. 

For all of these cases, "radically broken scale invariance" may prevail, in which the spectral 
index functions 

n,{k)-l = dlnV^^^Jk)/dlnk = 2{l + q-^) + q-^q'/{l + q)+ 2Cs, (2.1) 

q = dq/dlna = —qdq/dlji{Ha) , 
nt{k) = dlnVGw{k)/dlnk = 2{l + q-^) +2Ct, (2.2) 

can be relatively arbitrary functions of spatial wavenumber rather than constants. Here quantities 
such as q and q' are evaluated at the "time" Ha = k. The formula is motivated by the stochastic 
treatment of inflation fluctuations in the Hamilton Jacobi framework, in which quantum noise at the 
Hawking temperature H/2-k radiates from short distances across the decreasing {Ha)~^ boundary 
into a long wavelength background fleld. The post-inflation power spectra are parameterized by 

"-■■I- = ^iW^"'"^''""- ^OW=S^m2^f^'-. (2.3) 

Of course the utility of these expressions depends upon the correction factors C^^j, which are derived 
from the related Ug^t- Analytical forms for special cases can be derived, e.g., for uniform acceleration, 
and these show the Cg^t are typically small.* Much has been written on this subject. See, for 
example. Bond (1994), Lidsey et al. (1997), Wang et al. (1997) Schwarz et al. (2001), Kinney 
(2002), Habib et al. (2002), Martin and Schwarz (2003), Peiris et al. (2003), Leach and Liddle 
(2003), and references therein. The accurate path to the spectral indices is to take logarithmic 
derivatives of full numerical calculations to get the Cs,t- One can certainly invent cases in which 
the Cs^t are not small. However, provided q does not change too rapidly it is reasonable to use these 
formulas as guides. They show that tilt mostly depends upon how far the acceleration is below the 
critical value of unity. For qra —1, a substantial scalar tilt can come from the second term, yet no 
tensor tilt, as in natural inflation. 

Deviations from the power law model are not just expected, they are necessary, since q{k) 
must have passed from negative to positive to have created matter from the vacuum energy housed 
in the inflaton. The simplest form of braking of the acceleration is the running index, ns[k) = 

*The stochastic inflation technique uses "the H/(2ti) at fe = Ha WKB approximation", writes eq.JJSl &s a 
function of H , q and derivatives, and takes a logarithmic derivative wrt Ha in place of k. Analytical corrections 
invariably involve Hankel functions and their asymptotic expansions. No slow roll ((1 + g) ~ 0) restrictions are 
needed in these approaches. In the HJ formulation, H{(f)) and q{(j)) are treated as functions of the inflaton fleld, and 
satisfy the 'reduced Hamilton- Jacobi equation' relating H{(t>) to the potential V{4>): H^ — Hgji/{1 — (5+ 1)/3), where 



^i«=i7^and(l + g) = ^ 



ainH 



^ Mi. 



The extra piece in Us —nt is q ^q /(I + g) = — g ^^^ g2'^ 1 responsible 



for deviations in the two indices that can be signiflcant near (1 + g) ~ 0. 
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ns{kn) + [dns/dlnk{kn)] ln(/c/fc„). Such a form can be expressed in terms of coefficients g„, g^, q'^ 
in an expansion of q{k). 

A priori it may seem unlikely that a marked change in the expansion rate or acceleration would 
just happen to be in the 3-decade window of /c-space accessible to our CMB observations, since it 
maps into a relatively small patch on the inflaton potential surface because of the (l + (?) suppression 
factor in j^d(f> = ±^"-'^(1 + q)^'^dln{Ha). That is that q'^ would be small. However in (/)-space, 
this CMB window is not very far from (pend defining the acceleration/deceleration boundary, hence 
the q rise to zero must be reasonably rapid in </>. Even so, for most inflation models, the rapid 
change does indeed occur only near the end, suggesting special physics might have to be built 
in to accommodate large change earlier. Rapid acceleration changes, if present, would seem to 
be more likely a consequence of interaction with other field degrees of freedom rather than a 
result of inflaton self-interaction. Such hybrids involving two scalars interacting with either simple 
polynomial potentials (with second, third and fourth order terms), combinations of exponential 
potentials, and other simple forms, have long been used to show that constructing power spectra 
with mountains and valleys and even generating non-Gaussian fluctuations is possible in inflation. 

However even if more than one scalar field enters it is often possible to consider an effective single 
inflaton self-interacting through an effective single-inflaton potential over the observable scales. This 
is because the fields first settle into gorges on the potential surface, then follow the gorge downward 
towards the local minimum along a single field degree of freedom, (/>||, to be identified with the 
inflaton. The other degrees of freedom, (j)±, are 'isocurvature' degrees of freedom. Usually, the 
faces rising up from the gorge will be sufficiently steep that the inevitable quantum noise that 
excites motion up the walls quickly falls back, leaving no usable isocurvature imprint, effectively 
making those direections irrelevant (although curvature in the trough can lead to complications in 
the kinetic energy piece of the inflaton degree of freedom) . The single-inflaton expressions in terms 
of q would prevail. 

To get observable response often involves invoking an instability, with negative transverse com- 
ponents of the mass-squared matrix, d'^V/d<j)id(f)j, leading to an opening up of the gorge or its 
bifurcation. Tuning the location of such a structure to the window on the potential surface we can 
access may seem to be unpalatably precise. This is perhaps mitigated by relating it to a waterfall 
of sudden q change to trigger reheating, termed hybrid inflation. Control of residual defects left 
from the potential was always an issue, but such residual subdominant components are worthwhile 
to hunt for in the data. Although these multiple field models would have significant deceleration 
occurring in a waterfall phase, to have a spectrum with many sharp features littering the CMB 
range parameterizing a complex braking pattern seems very baroque indeed. (These multiple field 
situations are the ones where the simple spectral index formulas in terms of q are likely to have 
the largest corrections.) CMB phenomenologists should constrain such possibilities anyway. We 
shall see that the prognosis for constraining even such radical braking is reasonable with upcoming 
CMB experiments. 

The richness of infiation theory has expanded considerably with the emphasis on higher dimen- 
sions, brane-ology and stringy cosmology. So far to the extent there are predictions they fit within 
the basic infiation phenomenology as applied to CMB analysis of the sort described here. It would 
be nice if a smoking gun pointing to a uniquely stringy culprit will be found theoretically, and in 
the data. 
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III. Parameter Estimations, Now and Future 

A Parameters from the Jun03 Data 

Fig. El shows visually the one and two sigma constraints for the minimal inflation parameter set 
derived from the Jun03 CMB dataset using the Monte Carlo Markov Chain method. These results 
were also reported in Readhead et al. (2004) and are very similar to the Mar03 results given in 
BCP. BCP showed that the MCMC results were also in good agreement with those obtained using 
fixed C^-grids. Table ^ gives projected one-sigma MCMC errors. Priors applied to augment the 
CMB data with other information range from weak ones to stronger constraints from the expansion 
rate (HST-h), from cosmic acceleration from supernovae (SNl) and from galaxy clustering, grav- 
itational lensing and local cluster abundance (LSS). We show results in the table for CMB+weak 
and CMB+weak+LSS, with a flat i^tot = 1 prior also imposed. 

Fig.[7|and Table ^ address the level at which the WMAPl and Jun03 CMB data would prefer a 
running spectral index dng/dln k at about the two-sigma level. The projected distributions depend 
upon prior choices, e.g., Fig.EJshows that if one restricts how large tq is allowed to be, —dng/dln k 
would not be as large, accounting for the differences in the result here with those of the WMAPl 
team's analysis (Spergel et al. 2003). They highlighted how the non-CMB Lyman-alpha forest data 
in conjunction with the CMB suggested a running index. More effort is required to demonstrate 
that the forest estimates of power spectra are reliable enough to apply to this problem. 

Of course in statistics one should just ignore two sigma indications, and especially here when 
an extra cosmic parameter is added which has strong degeneracies with others in the basic minimal 
inflation mix. Fig. EJ shows the conspiracy driving the indications of running index from the Jun03 
data alone. The two low i anomalies and the slightly lower power at high i would prefer to 
bend the best-flt uniform acceleration model downward. By itself, WMAPl does this, and the 
addition of other experiments just takes this tendency and adds to it. However, BCP showed 
that the pre- WMAPl Jan03 compilation of the data that included DMR and Archeops also had a 
distribution that preferred negative [dns/dlnk{kn)], though with less statistical signiflcance than 
the post- WMAPl Mar03 set. 

Because theorists like to theorize about low significance results in anticipation they might even- 
tually emerge at high significance, much renewed discussion and many papers have now been written 
on whether the low i anomalies or the combination of low and high i anomalies indicate new physics. 
Within the context of inflation models, this involves arranging for q to change. If only low i is the 
target, a scale is built associated with a target k, e.g.. Bridle et al. (2003), Contaldi et al. (2003). 
Topology is another mechanism, building in a characteristic horizon-scale size to discretize fc-space, 
with just enough inflation to make the Universe just so big but no bigger. Trying to solve the high 
and low i anomalies with the same mechanism utilizes the running index, or would need to build 
in a mix of scales. 

Given this baroqueness, it is useful to explore the sensitivity of the running index distributions 
to cuts in the data. For example. Bridle et al. (2003) found [dns/dlnk{kn)] of — 0.04ib0.03 using all 
multipoles of WMAPl, the Jun02 versions of CBI and VSA and the Jan03 version of ACBAR, along 
with 2dF. When they excluded i < 5 from WMAPl, this dropped to -0.015 ± 0.03. To test this 
sensitivity further, we have marginalized over the i = 2,3,4,21,22,33 multipoles of the WMAPl 
data, which have "anomalous" bandpowers. We flnd —0.062 it 0.043 compared with —0.088 it 0.041 
with no such cuts. Just removing i = 22 gives —0.082 it 0.042. For this exercise, we did not 
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Figure 6: The state of parameter estimation using the Jun03 data compilation is illustrated by 
the one and two sigma contour regions when all but the two variables shown are marginalized. 
The scalar spectral index was not allowed to run. The outer contours are for WMAPl alone. 
The ID probability distributions for each single variable shown at right gives means and 1-sigma 
errors listed in Tabled This illustrates how the current data to higher £, predominantly driven by 
Boomerang, CBI and Acbar, sharpen the WMAPl results by breaking partial degeneracies. The 
priors applied were fi/; = and the weak /i-prior. In the figures, the ag-dominated LSS prior was 
included as well. The table shows the extent to which this sharpens parameter determinations. 
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Figure 7: Similar to Fig. El except the scalar spectral index is allowed to run. Note the ag and tc 
shifts indicated here and in Tabled 
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Table 1: Sample cosmological parameter values and their 1-sigma errors are shown, determined using MCMC methods after marginalizing 
over the other cosmological parameters and, for real data, the various experimental parameters. The first four columns show the state 
as of Jan04. WMAPl refers to the actual WMAP one year data, Jun03 is the compilation including the WMAPl data and recalibrated 
Boomerang, CBI, VSA and Acbar data, along with DASI and Maxima data. Means and standard deviations are given here. In all cases, the 
weak prior (0.45 < h < 0.9, age > 10 Gyr) is applied, and wq is fixed at —1, the cosmological constant case. The curvature parameter Clk 
is fixed at zero. The LSS prior agrees with current weak lensing data and agrees with most of the cluster determinations. A weak redshift 
survey constraint is also imposed. The 2dF "prior" uses the stronger 2dF redshift survey results, but not a ag prior. The primordial power 
spectrum index for scalar perturbations obeys ns{k) = ns{kn) + [dns/dlnk{kn)] (Ink/kn)- For the first set of numbers dns/dlnk{kn) is set 
to zero (no running of the spectral index). For the second set, it is allowed to vary. The last two columns show how the errors are expected 
to improve for WMAP with four years of data and Planck with one year of data, for the Fig. |^ simulations. 
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0.760 ±0.088 0.698 ±0.075 0.677 ± 0.050 0.669 ± 0.039 
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.02289 ± .00074 .02258 ± .00017 
.1113 ±.0070 .1165 ±.0014 
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Figure 8: The bottom panel translates the Q spectrum of Fig. ^ by mapping it from i into (per- 
pendicular) spatial wavenumber fe_L, where i ~ T^deck±, and also dividing out the target spectrum. 
It indicates the ratio of primordial observed Vc_{k) to target VY {k) in /c-space with a window 
function for each band. The window spillover into neighbouring bands is not large, so it shows 
where the anomalies lie for the specific best-fit ns=0.957 ACDM model. The downward bending 
curve is a shape that a running index of —0.09 would give. The low and high i downward-drive is 
evident. Breaking up the single low i bandpower into many bandpowers as shown in the middle 
panel highlights the lower I and £ ~ 20 "anomalies" . The right side of the middle panel shows ag , 
a broad-band power that probes cluster scales which reside near k~^ ~ 4h~^ Mpc, relative to the 
target model. The solid square is derived from the Jun03 data, with erg ss 0.85. The heavy open 
square to its left is the current best value for weak lensing (Hoekstra and van Waerbeke, 2004, 
private communication), which has evolved slightly downward from the Jun03 estimate shown as 
the leftmost light square. Many X-ray cluster estimates are lower (rightmost light butterfly). The 
heavy butterfly to its immediate left is the ag estimate from the SZ interpretation of the CBI, Acbar, 
BIMA high a. "anomaly" (e.g., Readhead et al. 2004). The top panel shows the recent SDSS data 
of Tegmark et al. (2004). Instead of dividing our Jun03 best-fit, which has some wiggles induced 
by the baryons, a best-fit wiggle-less "F-model" was used to highlight any obvious need, within the 
errors: the statistical answer is no (Tegmark et al. 2004). Galaxy biasing complicates the use of 
SDSS and 2dFRS, but there is no indication of ifiy running index within these LSS datasets. 



constrain the running index by any cosmological priors and allowed it to vary between ±0.3, which 
leads to a broadened tc distribution, as reflected in Table ^ 

Of course there should be a reason to justify removing or cleaning anomalies, e.g., that a source 
of systematic error is found or foregrounds and other residuals contaminate. In spite of much debate, 
there is no evidence that the low ^'s are low because of foreground contamination. However, when 
the quadrupole, octupole, etc., are obtained, the influence of the foregrounds should be reflected 
in the error bars. Slosar et al. (2004) improved the determination of errors on WMAPl at low 
I by marginalizing over foregrounds rather than using template subtractions. This leads to more 
power at low (. and better error determination. They improved the treatment of likelihood tails at 
low-£ over the "standard WMAPl prescription" of Verde et al. (2003), which we used to get the 
results given in our figures and tables. Both effects decrease the statistical desire for a downturn 
and Slosar et al. find that the ~ 2-0" effect drops to a ~ 1-a effect. 

Our conclusion, as in Bridle et al. (2003), BCP, Readhead et al. (2004) and Slosar et al. (2004), 
is that evidence for a running index in the CMB data is not compelling. To get the large values 
allowed by the data would require rather dramatic changes in the acceleration of the universe over 
what is actually quite a narrow range on the inflaton potential surface, manifested in a soft or even 
a radical breaking because of changing braking. 

B Forecasts of Parameter Precision 

The running index issue will probably be with us for quite a while, but forecasts are rosy for how 
well planned CMB experiments can answer this question: if there is a running index, it will be 
detected in the next generation of experiments, and if there is not it will be strongly constrained. 
This is illustrated in Fig. CUl Table El as weh as in Fig. Eland Tabled 

The precision for ground + WMAP4 would improve with larger sky coverage. For the numbers 
in the table, it was assumed that the primary spectrum beyond 2000 would be contaminated by 
secondary signals, but component separation should mitigate this, and the parameters are not 
sensitive to lesser cutting. Note how the larger baseline in (. significantly decreases the degeneracy 
and also the drift in the value of the running index one would estimate from the data relative to 
the target value. 

For the forecasts of Table [21 power spectra and their errors in Figs. |21 13 were calculated 
using "faster-like" algorithms of the sort we have appled to Boomerang. The idealizations make it 
"superfast" . The cosmic parameter errors were estimated using Fisher or curvature matrices about 
maximum likelihood values, rather than using Monte Carlo Markov Chains on the simulated data 
as in Tabled This means the parameters are treated as completely internal rather than completely 
external, as in MCMC. (See BCP for a discussion of the difference; the fixed grid approach was 
a hybrid, with amplitudes and experimental variables treated as internal, the rest external.) For 
both forecasting methods, the same Jun03 best-fit model with zero running index was used and the 
experimental parameters were essentially the same. The maximum likelihood drifts from the listed 
input cosmic parameter values of the target model depend upon the specific realization. Note that 
the errors of the Fisher-based forecast given here are quite similar to the MCMC values given in 
the last two columns of Table d 

Thus the superfast forecasts have been nicely validated by comparison with the MCMC results. 
They also give results compatible with what was actually obtained with Boomerang and WMAP 
when the real experimental specs were used. The forecasts can be made more sophisticated with 
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Figure 9: Forecast of one and two sigma contour regions for WMAP4 (black outer contours, light 
blue shading) and Planckl (red inner contours) show how the errors of Fig. [7] may improve in 
the future. Note some of the variables shown differ from those in Fig. [7] and there are changes in 
axes scales. Estimations were on simulated datasets generated using the no-running-index best- 
fit to the Jun03 CMB data. Only the 94 GHz channels were used for the 4 year simulation of 
WMAP4 and only the 150 GHz channels for the 1 year simulation of Planckl in these forecasts. 
No LSS priors were imposed. Table ^ gives means and errors for these two cases. The precision 
sharpens when all channels are brought to bear. Further, both Planck and WMAP are expected to 
observe for roughly double these periods, decreasing the noise component of the bandpower errors, 
with the sample-variance (cosmic-variance) component unaffected. One may therefore interpret 
these as conservative estimates, but the forecasts here do not include all of the extra complications 
associated with foreground separation. As Fig, ^J shows, anticipated ground-based experiments 
beyond the ones used in Fig. [7| will also have a powerful sharpening effect on precision. 
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some attempts at addition of foregrounds and residual signals and subsequent removal by param- 
eterized "power spectrum cleaning" . Tegmark et al. (2000) explored many aspects of the impact 
foreground contamination could have on forecasted cosmological parameter errors. Increasingly so- 
phisticated forecasts can still only be considered as partial steps towards the full mocks of a given 
experiment that one actually needs. 

Forecasting has a long history in the CMB, as a necessary ingredient for experimental proposals, 
and for showing feasibility of measuring new theoretical effects. Many realizations, experimental 
configurations and theoretical assumptions can be checked very quickly. This leads to an overwhelm- 
ing amount of information on the stages we expect to see between the data now and the Planck 
results, because there are so many polarization-sensitive ground-based experiments in various con- 
figurations either funded or proposed. Further there are many additional theoretical parameters to 
consider to further add to the information glut, yet turning them all on at the same time obscures 
what will happen in practice. 

In this paper, we have chosen to highlight only three forecasting cases, WMAP4, Planck and 
a fiducial high resolution ground-based experiment of modest sky coverage compared with what 
is possible. This reflects our experience that erring on the conservative side may reflect the real 
issue, which is how complex the analysis of the actual polarization data will turn out to be, and 
how much it will limit the precision we can obtain from the ground; and indeed from space. For 
polarization-targeting ground experiments, we can look forward to a wonderful set of developments. 
Listing some of the cases that have been considered gives an idea of the range. 

Apart from WMAP4 and Planckl with one channel, we have considered the following: WMAP2, 
WMAP4, WMAP8, using all 5 channels as well as the W channel with its 13' resolution adopted 
here, usually assuming 0.9 for f^^y - detector noise and beam sizes courtesy of Lyman Page; Planckl 
and Planck2, using either the 143 GHz channel with beam 7' alone, or with the 220 GHz, 5' and 100 
GHz, 9' polarization-sensitive bolometers (PSBs), or all together with the lower frequency HEMTs 
- from Planck "blue book" numbers, augmented by the most recent PSB numbers from Andrew 
Lange. With everything included forecasted errors do improve somewhat, but are often largely 
sample variance limited. 

We have considered forecasts for the Boomerang 2003 flight (the polarization analysis of the 
real data is heavily underway), and for Acbar (which is continuing to observe) - specs for current 
and subsequent observing seasons from John Ruhl. 

We have also forecasted for the South-Pole-based BICEP, using PSBs at 143 GHz, 40' and 100 
GHz, 60' observing 1000 sq deg in 260 days - numbers from Lange and Eric Hivon, with similar 
capabilities suggested for an experiment at another Antarctic site, Dome C, courtesy of Paolo de 
Bernardis. For QUaD (Quest mounted on DASI at the South Pole), we used PSBs at 143 GHz, 
4.0' and 100 GHz, 6.3' observing 200 sq deg in 260 days - specs from Lange and Hivon, Similar 
numbers are given for a separate Cardiff-based proposal. Prospects for BICEP and QUaD are very 
good for polarization and both are expected to be observing in 2005. See Bowden et al. (2004) for a 
full discussion of optimizing ground-based CMB polarization experiments, in terms of the tradeoff 
of sky coverage and sensitivity per pixel. 

The proposed QUIET experiment from Princeton, Chicago, JPL/Caltech, using new HEMT- 
MMIC array technology under development at JPL, would be mounted on the CBI platform in 
Chile, with large beams, 44 and 90 GHz, 42', over 2000 and 8000 sq deg, and small beams, 44 and 
90 GHz, 4', over 2000 and 8000 sq deg. Forecasts look very good for polarization and the first 
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phase could begin in late 2005. The South Pole Telescope numbers used only 220 GHz, 1.3', the SZ 
null channels - were from John Ruhl. We assumed a noise for polarization v2 times that in total 
anisotropy per pixel. For the Chile-based ACT, the beam is slightly larger than for the SPT. A 
deep mode of 100 sq deg was considered in addition to the 1000 sq deg we chose to highlight here - 
sensitivities from Lyman Page. In contrast to Fig. El this showed the power in the lensing-induced 
B-mode could in principle be detected. A fiducial CMBPol and an essentially cosmic variance 
limited all-sky survey at SPT/ACT resolution with very tiny noise have also been considered. 
Needless to say, for the latter the target parameters are recovered at the best-you-can-do level. 

We now turn back to the results shown in Table El The last set of rows show how the error bars 
open up when searching for more radical braking than the running index model gives. Consider 
the case when V,^{k) has a structure of unknown shape, as in radical broken scale invariance (BSI). 
For two given cosmological parameter sets, a Vc^{k) could be fashioned to morph one C\ into 
another. (See e.g., Souradeep et al. (1998), Wang et al. (1999) for a discussion of the role this 
degeneracy plays in parameter degradation.) Polarization information breaks this severe degeneracy 
because the peaks and dips of Q and C\ are in different locales. Parameterization is in terms 
of power amplitudes in a number of A:-bands of proscribed shape. 24 bands were chosen for the 
Table 121 case. Apart from the conventional banding in AP(^, we have expored the impact of band- 
colours (bands in An^), continuous wavelets, among others. The colour-banding makes more of 
a difference. However, although the errors determined are somewhat sensitive to the primordial 
spectrum band-type, band-placing and band-number, polarization does indeed nicely mitigate the 
effect of BSI-induced degeneracy for these planned experiments. Non-CMB information from LSS 
also helps to break degeneracies between cosmic parameters and 'P^(A;)-structure. Having tq from 
the low (. is important for breaking parameter degeneracies. 

As more parameters characterizing the inflation model are added, the precision continues to 
diminish unless near-degeneracies can be broken. For example, with a target value for Vqw IT^C, 
of 0.17, the Planckl realization of Fig. ^ shows a detection of the tensor B-mode is possible, 
which could lead to a good estimate of this amplitude. Getting the B-mode is very important for 
this. The specific example, determining 8 other cosmological parameters as well, gave a 0.135 it 
0.028 detection. With somewhat more optimistic noise forecasts, but allowing for the incomplete 
sky coverage mixing of E and B modes, Lewis (2003) finds that Planck should be able to detect 
primordial tensor modes at 95% confidence with greater than 95% probability if Vqw jT^C, ~ 0.03. 

We can conclude from exercises such as these on the experiments coming that parameters char- 
acterizing GW signals, mildly broken scale invariance associated with a running index, subdominant 
isocurvature components, and even radically broken scale invariance can be determined within the 
CMB data. 
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Figure 10: The forecasts of WMAP4 (green), Planckl (blue) and WMAP4 + ground-based 
ACT/SPT-like data (red) are shown compared with the target value (black dot). Projected one- 
sigma error bars are shown, and the one and two sigma ellipses which illustrate the correlations. 
The (magenta) cross shows the precision of the Jun03 data when the scalar spectral index is al- 
lowed to run, as it is for the simulations as well. The bigger (cyan) cross shows the state with 
the pre-WMAP Jan03 data. Large dedicated ground-based telescopes targeting high i with huge 
arrays of bolometers (e.g., ACT and SPT) or of HEMTs (QUIET), when combined with WMAP, 
should greatly increase parameter precision in the leadup to Planck. For this simulation, ACT/SPT 
experimental parameters were adopted, and the bolometers were assumed to be polarization sensi- 
tive to show their powerful impact on EE mode detection. The assumed coverage was 2.4% of the 
sky, 1000 square degrees. Increasing this would further improve the parameter estimates given in 
Table 121 since errors on many bands are sample- variance limited. Planck does so well because of its 
all-sky coverage, and will have a large impact on constraining multi-parameter deviations from the 
simple uniform-acceleration inflation models. 
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Table 2: Forecasted cosmological parameter values and their 1-sigma errors for the cases of Fig. 1101 
contrast what may be achievable for WMAP4, Planckl, and for future ground-based data (labeled 
ACT/SPT). The last column combines a polarization-sensitive bolometer-array ACT/SPT-like tele- 
scope (experimental choices as in Fig. Iiop with WMAP4. Increasing the assumed coverage beyond 
the adopted 2.4% of the sky would increase precision in these idealized forecasts. For WMAP4 and 
Planckl, the quadrupole was included and the useful sky coverage was chosen to be unity. When 
the quadrupole is excluded, the sky coverage is dropped to 90%, uJmu and T'gw/'P( are added to 
the parameter mix, and a tensor component and weak lensing effects on the power spectra are 
included in the target model, the errors grow modestly, usually within 20% of the sigmas listed. Of 
the 9 parameters, WMAP4 would determine 3 orthogonal combinations to ±0.01, 7 to ±0.1; with 
an SPT/ACT-like experiment, these rise to 4 and 8; and for Planckl, to 6 and 8. The lower rows 
illustrate the impact of including many more parameters to characterize ns{k) than the two in the 
running index model. Here 24 parameters defining the amplitudes of V(^{k) in 24 bands are added 
to the standard mix. The impact is not as severe as one might have expected since the polarization 
information breaks this severe degeneracy because the peaks and dips of C\ and C\ are in 
different locales than Q .In Souradeep et al. (1998), it was shown that parameter determination 
was significantly degraded if there was no polarization information when similar numbers of 'P^(A;) 
bands were added. 
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